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Abstract. We study the accelerated expansion of the Universe by using the kinematic 
approach. In this context, we parameterize the deceleration parameter, q(z), in a model 
independent way. Assuming three simple parameterizations we reconstruct q{z). We do the 
joint analysis with combination of latest cosmological data consisting of standard candles 
(Supernovae Union2 sample), standard ruler (CMB/BAO), cosmic clocks (age of passively 
evolving galaxies) and Hubble (H(z)) data. Our results support the accelerated expansion 
of the Universe. 
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1 Introduction 

Mapping the cosmic expansion history with accelerated expansion is one of the major chal- 
lenges faced by modern cosmology. Various observations, directly (Type la Supernovae, SNe 
la) and indirectly (CMB, BAO, lookback time) provide evidence for the recently observed 
accelerated expansion of the Universe [1-5]. The phenomenon of acceleration is usually at- 
tributed to some sort of dark energy [6]. The effect of dark energy is to change the sign 
of "deceleration parameter" q(z). Whether estimates of q(z) point towards accelerated or 
decelerated expansion, strongly depends on quality and quantity of the observational data 
at various redshifts. Therefore, one of the simple ways to understand the kinematics of the 
Universe is by phenomenologically parameterizing q(z). 

The phenomenological approach is advantageous since it does not rely on model specific 
assumptions like the composition of the Universe. It is assumed that our Universe is homo- 
geneous and isotropic at large scales and is described by a metric theory of gravity. This 
type of approach is explained in literature in a variety of ways [7—11]. 

Following the same line of thought we reconstruct the deceleration parameter using the 
Supernovae Union2 data set, CMB/BAO, Hubble parameter data (H(z)), and lookback time 
(LBT). Since the nature of the driving force of the Universe is still a mystery, therefore the 
choice of parameterization of q(z) is arbitrary. We use two two-parameter models of q(z) and 
and one one-parameter model to reconstruct the deceleration parameter. Riess et al. [12] 
showed using Gold (SNe la) data set that q(0) < at 99% confidence level, and there was a 
transition from recent acceleration to past deceleration. In this context the phenomenological 
approach to parameterize the q(z) is an easy way to find out q(0), and the transition redshift, 
zt, where expansion switches from being decelerated to accelerated. 

In addition to the data sets used in previous works, we have used the age of slowly 
evolving passive galaxies (lookback time) in our analysis. The addition of the LBT is signif- 
icant as it is complementary to the observations used in previous works and may help us to 
obtain more realistic constraints on the expansion history. 

The paper is organized as follows. In section II, we discuss the q{z) parameterizations, 
data and methodology. The results and discussion are explained in section III. 
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2 Model, Data and Methodology 



With the assumption that Universe is homogeneous and isotropic, and spatial flatness which is 
motivated by inflation and WMAP measurements [13], FRW metric describes the background 
geometry 

ds 2 = -dt 2 + a(t) 2 [dr 2 + rW] . 

Here a(t) is the scale factor. The cosmic scale factor is related to the redshift of free streaming 
photons in the usual way: a(t) = 1/(1 + z) . The expansion and deceleration rates can be 
defined as 
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where H is the Hubble parameter. We can write the Hubble parameter as 

1 + q(z') 



H = Hq exp 



1 + z' 



-dz' 



(2.1) 
(2.2) 

(2.3) 



Here Hq is the present value of the Hubble parameter. 

We parameterize the redshift dependence of q(z) in the following way, 

Qi(z) = qo + q\z , 

qil{z) = q 2 + q 3 
1 



qill(z) 
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The first parameterizations is a linear Taylor series expansion of the deceleration parameter 
around 2 = 0, where go and qi are its present value and its first derivative, respectively. The 
second parameterization has an advantage that it converges at high redshift, as expected. 
The last parameterization is such that it converges to the value ^ at high redshift, which is 
the value of the deceleration parameter at matter dominated epoch. 

In general the dark energy density parameter can be expressed as a function of redshift 

as 

1 + w(z') 



1 + z' 



Vt x {z) = fbcoexp 
where w(z) is the equation of state and is related to H(z) by 



dz 



w(z) 
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Substituting for H{z) we can express w(z) in terms of the deceleration parameter q(z) as[14] 

§(l + g (z))-l 
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2.1 Lookback time 



Most of the observations used for constraining cosmological parameters like SNe la, angular 
diameter distances etc. are distance based measurements. The LBT on the other hand is 
based on ages of distant galaxies. 

The LBT to an object at redshift z is defined as the difference between the present age 
of the Universe and its age at redshift z and can be calculated as 



f z dz 



)H(p) ' 



(2.7) 



where %(p) is the dimensionless Hubble parameter %(p) = H(p)/Hq and p are the parameters 
of the model. The observed lookback time t°£ s to an object at redshift Zi is defined as 
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where tff s is the observed age of the Universe. t(zi) is the age of the object defined as the 
difference between the age of the Universe at redshift Z{ and the age of the Universe when 
the object was born at redshift Zf. 



t{zi,p) = H Q 1 
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Unc = tQ bs — ti{zf) is the incubation time of the object. Since we don't know the forma- 
tion redshift of the objects in our sample, ti nc is treated as a nuisance parameter and we 
marginalize over it. 

From the definition of the LBT above, we can write 



t L (z,p) = H - 1 [ 
Jo 



dz 
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(2.11) 



Now, we can substitute for q(z) and write the lookback time for the three different 
parameterization. In the first case we get 



t L {z,p) = Hv l !\\ + z')^-^eM-1iz) dz 
Jo 

For the second parameterization the lookback time can be written as 



7(<?2 + 93 + 1, 53) - 7 ( 92 + 53 + 1, 



t L (z,p)=H - 1 q ^-'*- 1 

where 7 is the incomplete gamma function. For the last case, we get 

rz exp 
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To constrain the parameters we use the ages of 32 passively evolving galaxies in the 
redshift interval 0.117 < z < 1.845 [15]. We assume a 12% one standard deviation uncertainty 
on the age measurements. The value of the Hubble parameter Hq is fixed at Hq = 74.2 
km/sec/Mpc and the age of the Universe is taken to be t^ bs = 13.75 ± 0.13 Gyr. 

The likelihood function is defined as 

C oc exp 
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where \ 2 * s given by 

X 2 (Pi Hq, ti nc , t bs) 
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Here Ui is the uncertainty in the estimate of t(zi), <J t obs is the uncertainty in the estimate of 
to and tL(zi,p, Hq) and to(p,Ho) are the predicted ages. 

To marginalize over this nuisance parameter ti nc we define a modified log-likelihood 
function [16] 
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A = t L ( Zi ,p,H ) - [t$ s - t{ Zi )\. (2.19) 



We minimize the chi-squared with respect to the model parameters in the three cases to find 
the best fit values of the parameters. 

2.2 Hubble Parameter 

Measurements of the Hubble parameter as a function of redshift, H(z), can also be used to 
constrain the deceleration parameter. Stern et al (2010) [17], gave 11 measurements of H{z) 
in the redshift range 0.1 < z < 1.75. Further, Gaztanaga et al. (2009) [18], gave estimates of 
H{z) determined from line-of-sight BAO peak position observations. We use the combination 
of these two samples to constrain our parameters. We have a total of 13 data points [19-21]. 

For the three parameterizations the Hubble parameter can be written in terms of the 
model parameter as follows: 

H(z) = H exp(q lZ ) (1 + z ) 1+q °- qi 
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We find the constraints on the model parameters by minimizing the chi-squared function 
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2.3 Supernova Union2 data 

We use the Union2 compilation of 557 SNe la [22] for comparing the observed luminosity 
distance (derived from the distance modulus) with the theoretical luminosity distance. The 
distance modulus and the luminosity distance are related as 

t i th = m-M = 5log 10 - j ^ + 25 (2.24) 

Mpc 

Here m and M are the apparent and absolute magnitudes respectively. The luminosity 
distance is related to the Hubble parameter in a spatially flat Universe as 

r z dz 

d L (z) = (l + z) / — - (2.25) 
Jo H{z ) 

Substituting for the Hubble parameter in terms of the deceleration parameter we get 

d L = exp(gi) qf- qi 7 (<7i - *>,?!, + 1)) (2.26) 
-"0 

d L = exp(g 3 ) qs i92+q3) 7 U + 92, T ^-, q 3 ) (2.27) 
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in the case of the three parameterizations. Here 7 is the generalized incomplete gamma 
function. We find the constraints on the model parameters by minimizing the chi-squared 
function 

x Union2 {p, #0) - 2^ —2 ( 2 - 29 ) 



=1 a d L ,i 



2.4 CMB/BAO 



BAO refers to a length scale in the distribution of photons and baryons by the propagation 
of sound waves in the plasma of the early Universe and they can be treated as cosmological 
standard rulers. [5, 23]. The distilled parameter d z [24, 25] is defined as d z = > wnere 
r s (zd) is the comoving sound horizon size at the baryon drag epoch and Dy is the 'dilation 
scale' distance given by 

1 

\2n ^2 cz \ 3 



Dv(z) = 1(1 + z ) D a{zY^\ (2.30) 
where Da is the angular diameter distance given by 

, x 1 f z dz . 

d - (£) =(tW W) (2 ' 31) 

We can further use the measurement of the acoustic scale I a provided by CMB to define a 
ratio R z as 

1 - lAd '-(l + z.) D ^ T 'M (2.32) 
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The value of the ratio between the sound horizon at last scattering and at the baryon drag 
epoch is nearly 1.044 [25]. So we can write 



^^-(1 + ^) 1.044 ^1 (2 . 33 ) 

K z 7T JJy[Z) 

where z* is the redshift of recombination, r s (z*) is the size of the sound horizon at last 
scattering and Da(z*) is the physical angular diameter distance at the decoupling surface. 
We use six data points from SDSS LRG, 6dFGS and WiggleZ surveys. Three data points at 
high redshift (z = 0.44,0.6,0.73) are from WiggelZ survey [24]. Two data points at redshift 
(z = 0.2,0.35) are from SDSS LRG survey (Percival et al. [26]). The (z = 0.106) 6dFGS 
data point at low redshift is taken from Beutler et al. [27]. We derive the values of 1/R Z 
at 6 redshift points and find the corresponding errors by using the data and the correlation 
coefficients provided by Blake et al. [24]. We use ^=302. 09 ±0.76 [28]. 
The chi-squared function is given by 

X 2 CMB/BAO = (f° bS ~ f th f C ~ l (f° bS ~ f H ) ( 2 - 34 ) 

where f = and C is the covariance matrix evaluated using the correlation coefficients. 
The combined \ 2 is given by 

Xcombined ~ X Union2 + X Hubble + X LBT + X CMB j BAO (2.35) 

2.5 Estimating parameter errors 

The method for estimating errors on parameter values is given below. The \a marginalized 
uncertainties on the two parameters and their covariance can be calculated by using the 
following formulae [29]: 

2 _ ^ 2 _ a 2 _ C Ir. nc\ 
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Here L is the log-likelihood function, and qt and q{ are the best fit parameter values. 



3 Results and Discussion 



In this work we follow a model independent methodology to reconstruct the expansion history 
of the Universe. The significance of this kinematic approach lies in its simplicity because 
no dependence on the matter-energy contents of the Universe is assumed and further this 
approach does not demand any specific theory of gravity. This idea has been used in the past 
by many authors to prove the transition of the Universe from deceleration to acceleration 
phase. 

In 2002 by assuming a piecewise constant acceleration model with two distinct epochs, 
Turner and Riess showed that Universe will accelerate today if the transition redshift is 
fixed between 0.4 and 0.6 [9]. In a seminal work by Reiss et al. (2004), it was shown 
that Universe underwent transition from deceleration to acceleration by assuming the linear 
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parameterization of q(z) [12]. Further by using Gold SNe la data, Shapiro and Turner applied 
principal component analysis of q(z) and found very strong evidence (5<r) for the acceleration 
of the Universe in the recent past [7]. Elgaory & Multamuki (2006) [30], Gong & Wang (2006, 
2007) [31], Cunha & Lima (2008) [32], and Guimares, Cunha & Lima (2009) [33] again used 
SNe la data to map the kinematic expansion history of the Universe. 

Rapetti et al (2007) [34] developed new kinematical technique to study the expansion 
history of the Universe. They used the parameter space defined by the present value of 
deceleration parameter and the jerk parameter, j (dimensionless third derivative of scale 
factor w.r.t. cosmic time). By using SNe la data and X-ray gas mass fraction measurements 
they measured qo = —0.82 ± 0.14. Recently Lima, Holanda & Cunha (2009,2010) [35] used 
the Sunyaev- Zeldovich effect and X-ray surface brightness data to study the kinematical 
description of the expansion of the Universe. 

In this work we also followed the same line of thought and reconstructed q(z) by following 
purely kinematic approach. 

• Parameterization I: qi(z) = qo + qi z 

The la, 2a and 3a confidence contours in the go — Qi plane are shown for the com- 
bined data sets, both with LBT and without LBT data in Fig lb and Fig la respec- 
tively. The contours obtained with joint analysis are very tight as compared to the 
ones obtained with SNe and galaxy clusters. The best fit values of the model param- 
eters and the zt for different data sets are displayed in the table 1 below. The best 
fit values with la errors for (CMB/BAO+Hubble+Union2+LBT) data set and for 
(CMB/BAO+Hubble+Union2) data set are (q = -0.332 ± 0.018, q x = 0.146 ± 0.011, 
z t = 2.27 ± 0.12) and (q = -0.371 ± 0.023, qy = 0.154 ± 0.011, z t = 2.4 ± 0.22) respec- 
tively. The value of w(z) at present with la errors for (CMB/BAO+Hubble+Union2-|-LBT) 
data set is -0.792 ± 0.017. 



Table 1: Best fit values for first parameterization 



Data Set 


X 2 /d.o.f 


qo 


Qi 


?(0) 


Zt 


CMB/BAO 


0.538 


-0.491 


0.202 


-0.491 


2.43 


Union2 


0.954 


-0.524 


0.859 


-0.524 


0.61 


Hubble 


0.765 


-0.347 


0.449 


-0.347 


0.77 


LBT 


0.491 


-0.202 


0.07 


-0.202 


2.88 


CMB /BAO+Union2+Hubble 


0.971 


-0.371 


0.154 


-0.371 


2.4 


CMB/BAO+Union2+ Hubble + LBT 


1.009 


-0.332 


0.146 


-0.332 


2.27 



Fig. 3a and Fig. 3c show the evolution of deceleration parameter and equation of state 
respectively, with redshift z. There is a transition from accelerated to decelerated phase. 
The transition redshift in this model is z t = 2.4 (with SNe la + H(z) + CMB/BAO) 
and the zt is lower down to 2.2 when LBT data is included in the joint analysis of 
SNe la, H(z) and CMB/BAO. This zt is very high as compare to ACDM prediction 
(zt = 0.66). The joint analysis with LBT gives present value of deceleration parameter 
as qo = —0.332. 

• Parameterization II: qn{z) = qi + 53 (jt^) 

Confidence regions (la, 2a and 3a) in qi — q% parametric plane for SNe la, H(z), 
CMB/BAO and joint analysis are shown in Fig 2a and Fig. 2b. The best fit values with 
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(a) Without lookback time (b) With lookback time 



Figure 1: The la, 2a and 3a contours in q$ — q\ for first parameterization. The gray 
contours are for the combined chi-squared. Green (continuous line) contours correspond to 
CMB/BAO, red (Dash-dot) curves are for lookback time, orange (Dotted) lines are SNe la 
and blue (Dashed) curves are for Hubble data. 



la errors for the (CMB/BAO + Hubble + Union2 + LBT) data set and for (CMB/BAO 
+ Hubble + Union2) are (q 2 = -0.595 ± 0.073, q 3 = 1.278 ± 0.042, z t = 0.87 ± 0.20) 
and (g = -0.526 ±0.028 , q 3 = 1.205 ±0.046, z t = 0.77 ±0.09) respectively. The value 
of w(z) at present with la errors for (CMB/BAO±Hubble±Union2±LBT) data set is 
-1.042 ±0.07. 



Table 2: Best fit values for second parameterization 



Data Set 


X 2 /d-o.f 


02 


<?3 




Zt 


CMB/BAO 


0.556 


-0.775 


1.578 


-0.775 


0.96 


Union2 


0.956 


-0.552 


1.335 


-0.552 


0.70 


Hubble 


0.761 


-0.411 


1.014 


-0.411 


0.68 


LBT 


0.495 


-0.325 


0.58 


-0.325 


1.27 


CMB/BAO + Union2 + Hubble 


0.947 


-0.526 


1.205 


-0.526 


0.77 


CMB/BAO + Union2 + Hubble + LBT 


0.957 


-0.595 


1.278 


-0.595 


0.87 



Similarly the joint analysis is further done by the inclusion of LBT for this parameter- 
ization, shown in Fig. 2b. The joint analysis with LBT shifts the best fit parameter 
values towards the higher side. The constraints obtained on the parameter values by 
the joint analysis are very tight as compared to the constraints obtain from the SNe la 
data and the galaxy cluster data sets independently. 

Fig. 3b and Fig. 3d display the variation of deceleration parameter and equation of state 
respectively, w.r.t. the redshift. The reconstruction of the q{z) is done by the joint 
analysis of SNe la + CMB/BAO + LBT + H(z) data sets. The central line is drawn 
with the best fit values of the model parameters. The transition redshift in this case 
is, z t = 0.87. 
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12 12 

(a) Without lookback time (b) With lookback time 

Figure 2: The la, 2a and 3a contours in q2~Q3 plane with all the data sets in case of second 
parameterization. The gray contours are for the combined chi-squared. Green (continuous 
line) contours correspond to CMB/BAO, red (Dash-dot) curves are for lookback time, orange 
(Dotted) lines are SNe la and blue (Dashed) curves are for Hubble data. 



• Parameterization III: qjjj(z) = ^ + rff^p The best fit values with la errors for 
the (CMB/BAO + Hubble + Union2 + LBT) data set and for (CMB/BAO + Hubble 
+ Union2) are {q 4 = -1.296 ± 0.024, z t = 0.61 ± 0.01) and (q 4 = -1.162 ± 0.03, 
zt = 0.52 ± 0.02) respectively and hence the present value of deceleration parameter 
equal to qo = —0.796 (See Figs 4b and 4a). The value of w{z) at present with la errors 
for (CMB/BAO + Hubble + Union2 + LBT) data set is -1.234 ± 0.023. 

The best-fit evolution of deceleration parameter and equation of state with redshift is 
shown in Fig 4c and Fig 4d respectively. The la error bar in this curve is very tight 
as compare to the two parameter model of parameterization. The transition redshift 
in this case is zt = 0.61 which is in agreement with A CDM model with in la level. 
The results for single parameter parameterization for the rest of the data sets are 
summarized in table 3 below. 



Table 3: Best fit values for third parameterization 



Data Set 


X 2 /d-o.f 


<?4 


q(0) 




CMB/BAO 


0.589 


-1.264 


-0.764 


0.58 


Union2 


0.957 


-1.14 


-0.64 


0.51 


Hubble 


0.805 


-1.115 


-0.615 


0.49 


LBT 


0.567 


-1.47 


-0.97 


0.71 


CMB/BAO + Union2 + Hubble 


0.952 


-1.162 


-0.66 


0.52 


CMB/BAO + Union2 + Hubble + LBT 


1.028 


-1.296 


-0.796 


0.61 



The summary of the results are the following: 
• We have shown that inclusion of LBT in the joint analysis is very important in the 
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Figure 3: (a) and (b) correspond to variation of q(z) vs z for combined cm-squared for 
first and second parameterization respectively with la and 2a level, (c) and (d) correspond 
to variation of w(z) vs z for combined chi-squared for first and second parameterization 
respectively with la and 2a level. For plotting (c) and (d) the value of O m o is taken to be 
0.3. 



measurement of cosmological parameters. This addition of age of passively evolving 
galaxies as a cosmic clock is completely independent and competitive with standard 
candles and rulers. 

• For all the three parameterizations considered here, the data favors a transition from 
the deceleration phase to acceleration phase of the Universe. Also, at present our 
Universe is accelerating. This is also supported by the value of equation of state of 
dark energy at the present epoch (see Figs 3(c), 3(d) and 4(d)). 

• The transition redshift z% strongly depends upon the form of the parameterization of 
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(a) Without lookback time 



(b) With lookback time 




(c) (d) 

Figure 4: (a) and (b), Likelihood vs parameter plots with all the data sets in case of third 
parameterization. The filled gray plot are for the combined chi-squared. Green (continuous 
line) correspond to CMB/BAO, red (Dash-dot) curves are for lookback time, orange (Dotted) 
lines are SNe la and blue (Dashed) curves are for Hubble data, (c) q(z) vs z for combined 
chi-squared for third parameterization with la and 2a level (d) w(z) vs z for combined chi- 
squared for third parameterization with la and 2a level. For plotting (c) and (d) the value 
of Q m o is taken to be 0.3. 

q(z). The transition redshift from decelerated to accelerated expansion (q(z) = 0) 
increases with the inclusion of LBT in the analysis except for the parameterization I. 
Similarly the present value of deceleration parameter become more negative with the 
inclusion of LBT except for the parameterization I. But the change in the value of go 
is small. 

• The linear parameterization, qi, predicts the transition redshift zt > 2 which is not 
compatible with the ACDM model even at 3a level. So this model is not reliable at all. 

Using SNe la data only, Mortsell & Clarkson (2009) [36] showed that the expansion of 



Universe is accelerating at low redshifts even at > 12a level. In their recent work Xu et al. 
[37] put bound on model parameters of q(z) by using SNe la, BAO and observational H(z) 
data. Further Lu et al. (2011) [38] again tried to constrain the kinematic model by using the 
SNe la and H(z) data only. They show the two parameterized forms of q(z) clearly deviated 
from the ACDM model. 

It is natural to extend this work with addition of more observational data set like 
gravitational lensing as a new standard ruler, GRB's as a standard candles and age of globular 
clusters as a cosmic chronometers in the present work. 
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